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1.  INTRODUCTION 


Gas  tuibine  engine  vibration  is  influenced  by  mechanical  unbalances  of  various  engine  components 
in  addition  to  other  mechanical  faults  occurring  in  gear  boxes,  bearings,  piunps,  gears  per  se,  compressor 
rotors  and  stators,  as  well  as  shafts  (Zabriskie  1974;  Allwood  and  Christie  1991;  Kerfoot,  Hauck,  and  Palm 
1973).  The  vibration  sensitivity  to  each  foregoing  component  varies  widely  and  cannot  be  understood 
unless  each  individual  element  is  analyzed  separately  and  then  eventually  synthesized  coUectively.  It  is 
now  more  and  more  necessary  to  predict  accurately  both  the  local  and  global  performance  of  gas  turbine 
engines  for  diagnostic  purposes.  Although  in-depth  analysis  can  be  made  of  each  component  via  finite 
element  analysis,  these  methods  are  very  time  consuming  and  tedious  and  require  extensive  resources. 
A  quick  look  at  the  problem  from  a  lumped  parameter  perspective  is  sometimes  sufficient  to  understand 
the  diagnostic  situation  at  hand.  In  particular,  it  behooves  us  to  access  fimdamental  concepts  quickly  in 
order  to  appreciate  the  acceleratimis  and  forces  of  the  various  components  in  a  systematic  but  yet 
methodical  approach. 

Hence,  the  purpose  of  this  report  is  to  view  the  vibrational  phenomena  of  the  gas  turbine  engine  of 
an  Ml  tank  fiom  a  fimdamental  and  elementary  approach,  to  gain  familiarity  with  the  physical  components 
and  their  relation  to  each  other,  using  fiee-body  analysis  and  fimdamentals  of  engineering  concepts.  It 
is  with  certainty  that  a  unique  insight  can  be  gained  by  this  <q)pioach  that  can  never  be  achieved  via 
complex  numerical  techniques,  at  least  not  at  the  outset. 

It  is  quite  clear,  however,  that  in-depth  distributed  analysis  can  naturally  ensue  for  a  particular 
component  or  assembly  if  the  fimdamental  approach  reveals  that,  in  fact,  a  problem  exists  for  further 
investigation. 

2.  APPROACH 


In  appreciation  for  a  vibration  analysis  of  the  gas  turbine  engine,  we  should  first  establish  elementary 
principles  and  analogies  for  various  parts  of  the  overaU  structure.  This  can  only  be  done  through  an 
examination  of  first  principles  applied  to  elementary  example  components. 

Let  us  first  examine  the  mass-spring  system  as  follows: 
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Newton’s  Method 


Assuming  the  restoring  foree  is  in  the  elemental  spring,  we  have 


IFx  =  -kX, 


(1) 


such  that  the  force  is  acting  from  right  to  left. 


(1)  If  we  displace  the  mass,  m,  an  amount  X  positive  from  left  to  right,  then  X,  X,  F  arc 
also  assumed  positive. 


(2)  Using  Newton’s  Second  Law 


SF  *  mX, 


(2) 


=>  -kX  =  mX, 
or 


(3) 


X  +  JLx  »  0. 

m 


(4) 


The  solution  is  given  in  the  Appendix  as 


X(t)  =  X^cosmt  +  Z^sinwt, 
(0 


(5) 


where  ©  * 


_  1 
k 


m 


rad/s  (the  natural  circular  frequency). 


If  we  employ  the  energy  method,  then 
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(1)  The  energy  stored  in  the  spring  is 


(6) 


(2)  The  kinetic  energy  of  the  mass  is 


(7) 


(3)  Since  there  is  no  dissipation  of  energy, 

T  +  V  =  Constant 


ImX^  +  IkX^  =  C. 
2  2 


(8) 

(9) 


Hence,  if  we  differentiate  equation  (9),  we  have 


X  [mX  +  kx]  *  0. 


(10) 


Since  X  ^  0  for  all  time,  then  equation  (10)  reverts  to  equation  (9). 

Typical  of  gas  turbine  engines  are  rotating  shafts,  such  as  quill  shafts,  sun  gear  shafts,  and  shaft 
assemblies  transmitting  engine  torque  from  various  stage  rotors  to  gear  shafts.  Hence,  consider  a  disk  of 
mass  M  and  radius  R  at  the  end  of  a  "weightless"  shaft  fixed  at  the  other  end.  The  disk  vibrates  such  that 
the  restoring  torque  is  k^  <|)  in  which  the  torsional  spring  constant  becomes 


(11) 
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For  a  circular  cross-section,  shaft  of  radius  r. 

Ip  =  ^Tcr'*  (polar  moment  of  inertia)  (12) 

G  -  Shear  modulus  of  shaft  material. 

If  we  displace  the  disk  through  an  angle  <(),  the  restoring  torque  becomes  ^  ORjosite  in  direction  to  <|). 
Using  Newton’s  Second  Law, 


(13) 


such  that 


1  9 

Jq  =  _  MR'^  (mass  moment  of  the  disk  about  its  center  of  mass). 


(14) 


and 


^"^0  =  (torque  acting  on  disk  which  is  opposite  positive  direction  of  <|>).  (15) 


We  obtain, 

-k^(|)  -  Jq^, 
or 

+  k^<t>  =  0  =  lMR2<t;  +  ^<1)  =  0, 

2  L 

or 


(16) 


(17) 


$  + 


TcGr'*  , 
- o 

MR^L 


0. 


(18) 


The  solution  of  Equation  (18)  is  given  in  the  Appendix  as 
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where 


(j)  =  (l>o  COSCDt  + 


_ sincDt, 

Q) 


(0  = 


.12:1]^  rad/s. 
MR^L 


(19) 


(20) 


Considering  other  turbine  engine  components  analogizing  a  pendulum  such  as  the  connection  of  the  inlet 
housing  assembly  to  the  air  diffuser  housing,  we  have, 

ZMq  =  -mgLsinO,  (21) 

displacing  mass  counterclockwise  as  positive  such  that  6,  6,  and  M  are  positive  counterclockwise. 

Using  Newton’s  Second  Law, 

IMq  =  IflO  *  -mgLsinO,  (22) 

where 

lo  *  mL^,  (23) 

(moment  of  inertia  of  its  mass  with  respect  to  axis  counterclockwise  normal  to  the  view).  Thus, 


or  for  0  small. 


Therefore, 


mL^0  +  mgL  sin0  =  0, 


sin0  »  0. 


0  +  1.0  =  0. 

L 


(24) 


(25) 
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Equation  (25)  was  solved  in  the  Appendix.  The  solution  is  given  by, 


0  =  0n  cos  ©t  + _ sin  ©t, 

© 


(26) 


where 


©  = 


.  .  1 

_g 

L 


rad/s  (natural  circular  frequency). 


(27) 


If  we  apply  the  energy  method  to  the  same  component,  then  let  us  assume  that  the  change  in  potential 
energy  in  the  system  is  due  to  mass  m  moving  upward  (gaining  energy)  in  a  gravitational  field. 

V  =  mg(L  -  L  cos  0).  (28) 

The  instantaneous  velocity  of  the  mass  is 

V  =  L0.  (29) 

Thus  the  kinetic  energy  is 

T  =  ymv^  =  .im(L0)^  (30) 


Since  there  is  assumed  no  loss  in  energy  (no  damping), 

(31) 

(32) 


T  +  V  *  constant. 


or 


.lmL^0^  +  mgL(l  -  cos0)  *  C. 


If  we  differentiate  this  expression  with  respect  to  time,  we  obtain 
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(33) 


mL^60  +  mgL6  sin0  =  0, 


or 


0  [mL^O  +  mgL  sin  0]  =  0. 


(34) 


Since  0  0  for  all  time  (no  motion),  then 


mL^0  +  mgL  sin0  =  0, 


or  for  0  small. 


0  +  i.0  =  0, 

L 


(35) 


(36) 


whose  solution  is  given  as  for  Equation  (25). 

To  analyze  the  inlet  guide  vane  (IGV)  slotted  lever  associated  with  the  IGV  system,  we  can  anaingj/f 
the  situation  through  a  mass  spring  system  by  the  use  of  Rayleigh’s  Method  to  include  the  effect  of  the 
mass  of  the  spring. 

If  we  let, 

Z  -  which  locates  a  point  on  the  spring 
u  -  which  is  the  definition  of  the  spring  at  Z, 

then. 

Kinetic  energy  =  T  = 

Stored  energy  =  V  = 

we  assiune  u  = 


>  J[<(pdz)0' 


IkX^ 

2 


ix 


(37) 

(38) 

(39) 
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then, 


T  =  ImX^  + 

2  2 


M 


m 


X^. 


Also,  since  m  =  pt,  it  follows  that  (O  = 


{2  2  2 

3  J 

nl 

K 


xf 

M  + _ 


,  the  Rayleigh  frequency. 


(40) 

(41) 


(42) 


Thus,  one-third  of  the  mass  of  the  spring  is  added  to  the  mass  M  to  account  for  the  effect  of  the  spring 
mass. 


Furthermore,  if  we  could  extend  the  IGV  slotted  lever  analogy  to  include  a  forced-simple  mass  spring 
system  where  sin  Qt  is  the  force  and  k  is  the  spring  constant,  then 


ZFx  =  m  X , 


(43) 


XFx  =  Fq  sinflt  -  kX, 


(44) 


so 


mX  +  kX  =  Fq  sincot. 


(45) 


or 


x  +  sinQT  =  ^  Ji  sinflt, 

mm  k  © 


(46) 


or 


X  =  5^  (0^  sin  fit. 


(47) 
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where 


(48) 


or 


X  +  sinQt. 


(49) 


The  steady-state  solution  as  given  in  the  Appendix  is 


X,  =  5, _ t _ sin^2t. 

S  S  o 

1  - 


(50) 


In  the  system  shown  in  Figure  1, 


a  10  in,  b  —  15  in,  mg  —  5.0  lb,  and  k  —  50  lb-in.  We  can  then  determine  (a)  the  natural  frequency,  and 
(b)  the  steady-state  solution,  i.e.. 


•  S  (Torque)o  such  that  Jq  =  mb^ 

(51) 

=>  mb^e  +  ka^e  =  Fjbsinflt 

(52) 

..  F,b 

0 - -  6  = - sinQt 

mb  2  mb2 

(53) 
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Noting  that 


0)2  =  Jii!.  =  50  X  (10)2 

_L  .  (15)2 
386 


1717.33 


(54) 


where,  g  ~  386  in/s^ 


{  Note: 


m 


-  [386  in/s 


(55) 


to  =  41.44  rad/s;  Fj  =  25  lb;  =  50  rad/s. 


(56) 


Thus, 


AMPL  =  0(,  =  02 . 


F,b 


25  X  15 


1  - 


ka" 


1- 


50(10)' 


1  - 


50 


41.44 


(57) 


or 

00  =  0.164  radian  =  9.428  degrees . 

Note:  r  ^  frequency  ^ 

^  Free  vibration  frequency^ 

To  analyze  the  torque  distributed  to  the  quill  shaft  in  the  Accessory  Geaitx)x  Module,  let  us  consider 
a  2-in-thick  steel  disc  having  a  radius  of  8  inches  subjected  to  an  oscillatory  torque  of  1,000  sin 
lOOrt  t  in/Ib.  The  steel  shaft  fiom  which  the  disc  is  suspended  is  fixed  at  its  upper  end,  has  a  2-in 
diameter,  and  may  be  considered  weightless.  We  can  easily  determine: 

(1)  The  natural  frequency  of  the  system, 

(2)  The  steady-state  angular  amplitude  of  motion  of  the  disc, 

(3)  The  maximum  oscillatory  shear  stress  in  the  shaft. 
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Thus, 


Mass  of  Disk  =  pV  = 


0.3 

386.4 


X  jc(8)^  X  2. 


(58) 


Mass  Moment  of  Inertia  =  0.3122  %  s^  •  in,  i.e.. 


J  =  1mr2  =  Ix  0.3122(8)2  =  9.99  lb.  in  s^. 
2  2 


(59) 


Shaft  Spring  Constant: 


Ke 


12  X  10®  X  1.57 

■  -  ss 
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1.256  X  10®  in/lb/rad. 


(60) 


where 


Polar  Moment  of  Inertia  is 


1  = 

P  ~W  ~  “32”  “ 


(61) 


since 


to 


2 


1.256  X  10® 
9.99 


=  0.126  X  10® 


(62) 


=>  CD  =  354  rad/s  or  f  =  =  .Hd.  =  56  cps  (Hz). 

2k  2n 


(63) 


Also, 


where 


(64) 


=  lOOn  »  314.16  rad/s 


(65) 
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and 


r  =  =  0.887 

(o  354 


(66) 


1,000 


1.256  X  10^  1  -  (0.887)^ 


=  3,734  X  10' 


(67) 


Thus, 


Shear  stress  =  x  = 


TflC 


K%)  X  1 

1.57 


Ib/in^, 


(68) 


where  Tq  =  torque,  and 

C  =  distance  from  neutral  axis  to  extreme  fiber. 


Hence, 


1.256  X  10^  X  3,734  x  10"^  x  1 

Tst 


2,987  Ib/in^. 


(69) 


The  situation  of  the  tie  rods  fastening  the  forward  assembly  to  the  rear  assembly  can  be  typified  as 
a  cantilever  beam  with  a  200-lb  weight  at  its  tip.  An  exciting  force  equal  to  F  =  5(X)  sin  20  7t  t  acts  on 
the  mass  along  its  vertical  centerline.  We  can  then  determine 


(1)  The  natural  frequency  of  the  system, 

(2)  The  steady-state  amplitude  of  the  motion  of  the  mass,  and 

(3)  The  maximum  bending  stress. 


Hence,  we  have. 


M  = 

386.4 


=  0.517  lb„ 


W 

T 


(70) 
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and 


with 


and 


^  3EI  3x30x10^x  0.667  ,,,. 

K  =  __  = - ^ - Ib/in 

(24)3 


I  =  ^bh3  =  _i_  X  1  X  (2)3  =  0.667  in^ 
12  12 


(71) 


(72) 


K  =  4,342  Ib/in . 


(73) 


Therefore, 


<{?■  = 


_K 

M 


4,342 

0.517 


=  8,388 


(74) 


Thus, 


where 


=>  CO  =  91.59  rad/s,  f  =  .^  =  14.5  Hz . 

2ic 


500 

4,342 


1  - 


20  Jt 


91.59 


0.217  in. 


(75) 


(76) 


and  the  steady-state  displacement  =  0.217  sin  20  n  t. 


If  the  bending  stress  at  the  support  =  o  =  , 


(77) 


=>M  =  (kx)l  =  (4,342  X  0.2 1T)  x  24  =  22,613  in/lb  (78) 
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whereas 


22,613  X  1 
0667 


=  33,900  Ib/in^  = 


TnC 


(79) 


The  IGV  operation  occurs  when  fiiel  from  the  electromechanical  fuel  system  enters  the  IGV  actuator 
and  causes  the  piston  to  move.  The  accompanied  movement  might  be  typified  by  a  mass-spring-damper 
system  such  that 

Exciting  Force  -  Fq  cos  £2t 
Fq  -  Amplitude 

-  Forcing  Frequency  -  rad/s 

such  that 


MX  +  cX  +  kX  =  FQCosfit, 


(80) 


or 


X 


—  cosGt. 
M 


(81) 


The  transient  solution  is  (see  the  Appendix) 

Xj,  =  cos  (©gt  -  ((>)• 


(82) 


This  response  approaches  zero  after  a  brief  time,  therefore,  only  the  steady-state  response  exists,  i.e.. 


where 


Xp  =  Xg  cos(f2t  -  <t)i). 


Xg  -  dynamic  amplitude  = 


<1  -  -  (25r)- 


(83) 

(84) 
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and 


_  Forcing  frequency 
CO  Undamped  natural  frequency 


(85) 


and 


E  s  ^  in  this  case  E 

(D 


c 

2m© 


(86) 


also, 


Fa  /  \ 

6j  = _  ^Static  displacement  of  M  due  to  a  force  Fq  j. 


(87) 


Hence, 

(1)1  =  tan"^_^ 
1  -  r 


(88) 


which  is  the  phase  angle  between  the  force  and  displacement. 


At  resonance  r  =  1, 


Xo  1 


90° 


As  r  ->  oo. 


<|)i  ^  139°. 


(89) 

(90) 


Identifying  the  dynamic  force  amplitude  transmitted  into  the  support  (Figure  2), 
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we  have 


Fj  =  CX  +  kX  =  Fjq  cos  (Qt  -  <t)|  +  <t)2). 


(91) 


Since 


6i  cos  (tit  -  <t>i) 

[(1  -  r2)2  +  (2^r)2j 


^TO  “  ^0 


1  +  (2^)2 


L(1  -  r2)2  +  (2^)^ 


(92) 


(93) 


hence, 


Ftq 

T  =  The  transmissibility  =  = 

1  +  (2^rf 

Fo 

,(l  -  *  (2^rf . 

(94) 


"Isolation”  of  various  components  for  analysis  within  assemblies  is  common  for  all  gas  turbine 
engines.  If  we  consider  in  a  generic  sense  any  of  the  turbine  engine  modules  which  are  modeled  as  a 
spring-mass-damper  system,  then  we  have 


Xi  =  XjQ  cos  £2t  =  displacement  of  support  or  "box" 
X2  =  Absolute  displacement  of  mass  to  be  isolated. 

Equation  of  motion  of  mass  m  is 


X2  +  i.  X2  +  JL  X2  =  JL  Xj  +  ±  Xj  =  Xjo 

m  ^  m  ^  m  ^  m  ^ 


(D^ 


1 

+  (2^r)^]^cos(f2t  +  <[)), 


(95) 


where  the  solution  is 


^2  ~  ^20  (^1  -  <j)  +  ^2)  ®  XjQ 

1  .  (2^r)2 

,(l  -  r^r  *  (25r)^J 

cos  (qi  -  (|)  +  ^2)- 


(96) 
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The  ratio  of  the  amplitude  of  the  isolated  mass  m  to  the  amplitude  of  the  support  or  "box"  is 


1  -  (2^r)2 


1 

T 


,(l  -  rf  . 


such  that  this  ratio  plots  exactly  as  a  "transmissibility"  curve. 


(97) 


Considering  the  rotating  system  assembly  and  its  associated  platfonn,  if  the  assembly  and  platform 
are  displaced  and  released,  then  the  displacement  vs.  time  might  look  like  the  following  (Figure  3) 
(assuming  a  total  weight  of  80  lb). 


Figure  3.  Displacement  vs.  time  for  rotating  system  assembly. 

What  is  the  damped  natural  circular  frequency?  What  is  the  damping  ratio  ^7  What  is  the  undamped 
natural  circular  frequency? 

Let  us  proceed  as  follows: 


5 


=  In 


X 


n  +  1 


=  In  i£  +  0.5877  =  271^ 
1.0  ^ 


(98) 


so  that  ^  =  0.0935 . 
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From  X  vs.  t,  one  cycle  takes  0.1  s 


(Oj  ({OR 

/.  ®  ^ - -  =  62.8  rad/s.  (102) 

[l  -  ^2]^  [l  -  (0.0935)^)^ 

If  an  exciting  force  10  sin  30  Ji  t  acts  on  the  assembly  and  platform,  then,  what  is  the  steady-state 
displacement  and  what  force  is  transmitted  into  the  support  as  shown  in  Figure  4? 
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Since 


and 


and 


X  = 


.(l  -  rf  *  i2irf. 


1 

T 


(103) 


Fq  «  10  lb 


(104) 


since  (O^ 


k 

M 


then  k  =  Meo^  =  J2.  x  (62.8)^ 
386 


(105) 


=>  k  =  817.37  Ib/in 


(106) 


and 


5i 


10 

817.37 


0.0122  in;  r  =  ii  = 

(D  20ic 


1.5, 


(107) 


hence, 


X 


0.0122 


1 


+  (2  X  0.0935  X  1.5)^ 


T 


0.0095  in. 


(108) 


Also,  the  transmitted  force  is 


Ft'Fo 


1  +(2^r)' 


1 

IT 


[(l-r2)\(2^rp. 


6.5  lb. 


(109) 


If  we  look  again  at  the  IGV  operation  returning  us  to  the  mass-string-damper  system  as  in  Figure  2, 
but  with  a  rotating  unbalanced  external  force  due  to  perhaps  wearAear,  we  have 
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MX  +  CX  +  kX  =  me£2^cosQt 
(Note  that  the  horizontal  component  =  me  cos  Q  t.) 


(110) 


or 


X  +  -£X  +  —X  =  rncOt 

M  M  M 


(111) 


m  -  rotating  mass;  e  -  eccentricity  -  inches 
Q  -  angular  velocity  -  rad/s 

The  steady-state  solution  is 


where 


Xp  =  .^e 

P  M 


COS 


.(l  -  rf  *  {2Kr)\ 


(sit  -  #1), 


(112) 


r  = 


SI 
^  ’ 


(113) 


c 

2M©  ’ 


1  - 


(114) 


Taking  into  account  Force  Transmissibility,  we  have 


And  noting  that 


Ft-  -  kX  +  CX. 


r 


2 


.(>  -  rf  .  (25 


cos{flt  -  (j)), 


(115) 
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1 


(116) 


=»  X  = 


Ma)2 


,(l-rf  .(25rp. 


1  1 
7 


cos  (£2t  -  (j)). 


Now,  since 


Hy  —  l*TO  COs(Ot  “  (|>) 


1 

7 

(117) 


Also, 


T  R  =  ^*1*0 


1  +  (2^r)^ 


L(l  -  rf  .  {2%t)^\ 


(118) 


Fjo  =  meiJ 


1  +  {2%Tf 


(l  -  rf  .  (25r)^, 


As  a  specific  jqqilication  of  the  unbalanced  force  phenomena,  consider  a  typical  gas  turbine  engine 
rotating  gear  shaft  weighing  40  lb  with  an  unbalanced  torque  ^lied  of  0.5  in/lb.  Some  tests  indicate  a 
natural  frequency  of  1,000  ipm  and  a  damping  ratio  of  0.2,  Let  us  determine  the  steady-state  amplitude 
when  operating  at  1,200  ipm  and  the  amount  of  imbalance  force  transmitted  into  the  support. 

Letting  r  =  £  =  ^  =  1.2,  r^  =  1.44,  ^  =  0.2  2^r  =  0.48  (119) 


^  Q 


1,200  X  =  125.66  rad/s 
60 


(120) 


and 


me  =  —  =  0.00129  in/lb, 

386.4  in/s  2 


(121) 


where 

g  =  386.4  in/s  2  =  32.2  ft/s^ 


(122) 
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such  that 


40 

M  =  . . ^  =  0.1035  lb. 

386.4 


(123) 


Steady-state  displacement  gives  us 


X  = 


me 

TT 


0.00129 


1.44 


(l  -  rf  .  (24r)- 


0.1035  1 

[(l  -  1.44)^  +  (0.48)^]  ^ 


(124) 


X  =  0.0276  in. 


(125) 


Since  Ep  =  Force  transmitted, 


=>  Ft  =  metf 


1  +  (2^r)' 


(■  -  rf  .  {l%rY 


=  34.7  lb. 


(126) 


If  we  have  a  periodic  but  nonsinusoidal  excitation  of  an  engine  component  analogous  to  our  spring- 
mass  system,  then. 


where 


F(t)  =  Ao  +  'L  sins  £2t  +  Z  cos  flt, 

Sel  S=1 


271 

^  F(t)  sinsfltdt. 


(127) 


(128) 


2n 

^  F(t)cossf2tdt. 


(129) 
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For  the  previous  case. 


Jt  27t 

n 


(130) 


A.-a  F 


2n 


(131) 


B.  =  _^  2F/)  cos  s  Q  t  dt  + 

®  2;tJo  0 


2k 

TT 


0  •  dt  =  0 


(132) 


F{t)  -  Fo  + 


~  2Fo 

I>  — -  sin  s  £2  t  . 

1,  3,  5.  7  Jts 


(133) 


Then, 


X  +  (D^X 


2F, 


^  sin  s  Q  t . 


(134) 


whose  solution  is 


X  = 


sin  s  Q  t. 


(135) 


where 


ifl,  and  (0^  =  JL. 
(0  m 


(136) 


Most  of  the  vibrational  problems  associated  with  gas  turbine  engines  involve  torsional  analyses  of 
components,  viz.,  rotating  gear  shafts,  rotating  turbine  blades,  rotors,  nozzles,  bearings,  spacers,  shims,  etc. 
In  many  instances,  two-degree-of-freedom  studies  are  required  for  consideration  of  vibrational  diagnostics. 
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Lx)oking  at  a  typical  disk-shaft  configuration,  let  us  consider  the  following  generic  situation  of  two  disks 
in  torsion  on  a  shaft: 


Z  TorquCfl  =  J0  -kjGj  -  ^2  (^i  -  62)  =  ®i- 

Hence, 


Also, 


Jj  61  +  |kj  +  kjj  Gj  +  k2  Gj  =  0 . 


(137) 


ZT  =  J2  02  k2  [Gj  -  G2)  =  ^2^2- 


Hence, 


J2G2  +  k2G2  -  k2Gi  =  0. 


(138) 


Assuming 


Gj  =  Aj  cos  (cat  +  \ff) 
02  =  A2  cos  (cot  +  \|r). 


We  obtain  from  (137)  and  (138) 


‘•2  2 

—  - 


(139) 


For  kj  -  k2  -  k  and  Jj  =  2J2,  then  the  |Det|  =0  yielding  the  frequency  equation. 
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(140) 


or 


where 


solving 


(2k 


(141) 


(142) 


©1 


(oj  =  0.293  0)J 


(143) 


CDr 


=  1.707  0)0. 


(144) 


The  mode  shapes  may  be  obtained  from 


or 


-oa^Ai  + 


=  0. 


(145) 


(146) 


For  0)2  =  ©j  =  0.293  ©q 


All  (1  -  0-293) 


©r 


=  0.707 , 


'21 


0)r 


For  ©2  =  ©2  =  1.707  Wq 


Ai2  (1  -  1.707) 


(On 


=  -0.707 . 


^22  (fit 
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Consider  a  two-degree-of-freedom  system  having  a  forcing  function  Fj  sin  Q  t  acting  on  mass  mj. 
The  coupled  equations  may  be  written  in  the  matrix  form: 


'mi 

0  ■ 

Xi 

'  + 

1^11 

^12 

. 

Xi 

0 

m2 

X2 

^  J 

1^22  _ 

Fj  sin  Q  t 
0 


(147) 


For  a  steady-state  solution,  assume 


Xj  =  Xj  sin  SI  t  and  Xj  »  Xj  sin  t . 


(148) 


The  two  equations  above  become: 


-  ki2 

*^12  (^22  - 


(149) 


Using  Cramer’s  Rule, 


Fl 


k 


12 


(150) 


such  that  D  =  determinant  of  coefficient  =  (k^  -  -  kjj  and 


X2  = 


-Fiki2 

D 


(151) 
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If  we  let  cOj  and  co^  be  the  two  natural  frequencies,  then  D  may  be  written  as 


D  =  m2m2(a)i  • 


(152) 


where  (0]  and  002  are  found  by  setting  D  =  0. 

An  alternative  analysis  or  formulation  using  transfer  matrices  can  be  employed  for  obtaining  relations 
that  govern  the  motion  of  a  discrete  system  composed  of  liunped  masses  cmmected  by  massless  elastic 
parts.  The  properties  and  conditions  are  expressed  by  state  vectors  at  sections  or  points  immediately 
adjacent  to  the  sides  of  a  discrete  mass.  Specifically,  a  state  vector  is  a  column  matrix  which  contains 
the  components  of  the  displacements,  forces,  and  moments  at  a  point  or  station  adjacent  to  a  mass.  Such 
a  state  vector  can  then  be  transferred  to  another  location  by  a  transfer  matrix,  there  being  two  types.  A 
point  transfer  matrix  transfers  a  state  vector  from  a  location  on  one  side  of  a  mass  to  the  odier  side,  at  the 
same  designated  station  and  thus  is  a  transfer  at  a  point  A  field  transfer  matrix  transfers  a  state  vector 
across  a  spatial  distance  or  field  of  the  system  from  a  station  at  one  mass  to  a  station  at  another  mass. 

The  relations  resulting  from  the  use  of  transfer  matrices  lead  to  a  solution  in  which  the  natural 
frequencies  and  mode  shapes  may  be  determined  from  the  characteristic  equation  if  the  number  of  degrees 
of  freedom  is  small  or  otherwise  by  a  numerical  procedure  such  as  Holzer’s  method. 

Consider  a  general  mass-spring  system  that  is  restricted  to  move  in  the  horizontal  direction  only.  The 
letters  "L”  arxl  "R"  are  used  to  denote  the  left-hand  and  right-hand  sides,  respectively,  of  a  mass  station. 

For  spring  kj,  the  following  relation  can  be  written: 


(153) 


F  ^  ®  F 
j  j  -  1- 


(154) 


27 


In  matrix  form,  these  would  be  written  as 


Or,  more  concisely, 


where 


(155) 


(156) 


(157) 


are  state  vectors  for  the  displacements  and  internal  forces  at  stations  j  and  j  -  1,  respectively.  The  scalar 
matrix, 

(F]i  = 

is  the  field  transfer  matrix  which  relates  {v}|'  to  {v}j^_  j 


1  1 
k 


0  1 


ji 


(158) 


For  mass  mj,  the  following  equations  apply: 


(159) 


mjXj. 


(160) 
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Now  for  hannonic  motion  of  mj,  Xj  *  Aj  sin  (tot  +  (t>)  and  Xj=  Aj  sin  (©t  +  <]>),  so  that 


Xj  =  -CD^ 


(161) 


Equation  (163)  then  becomes 


Fj'  .  -mjrfxj*-  *  Fj*-. 


(162) 


In  matrix  form,  Equations  (162)  and  (165)  would  be  expressed  as 


or 


{v}f  =  [P],{v}{- 


(163) 


(164) 


Here  {v}j^  and  {v}|' 


are  state  vectors  for  the  displacements  and  internal  forces  to  the  right  and  left 


of  mass  nij,  respectively,  and 


0 


(165) 


is  the  scalar  point  transfer  matrix,  which  relates  {v}j^  and  {vlj'.  Substituting  Equation  (159)  into 
Equation  (164)  gives 


{v}f=[P]j[Fl3{v}f.,=[Qj.{v}^.,, 


(166) 
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where  (see  Equations  [158]  and  [165]) 


[Q]j  =[P]j[F]j. 


(167) 


1  0 

1 

f 

¥ 

-mto^  1 

j 

0 

iJj 

(168) 


(l  X  1  +  0  X  O) 
(-mo)^  X  1  +  1  X  o) 


I  k 

2  1 
-mo)  X  _ 


0 


X  1 

7 


'  ^  'Ji 


1 

f 

-mo)^  1 
V 


1 

k 


ni(o^ 


•  (169) 


This  procedure  can  be  continued.  Thus 


WMq],M5. 

where 

(<5]n-[PL[Fl„[P]„-,[Fl„_,-(p],[F],. 


(170) 


(171) 


Equation  (170)  expresses  the  state  vector  {v}^  at  the  nth  station  in  terms  of  the  state  vector  {v}q  at  the 
initial  station. 


Noting  that  Equation  (171)  represents  two  algebraic  equations,  and  that  usually  a  boundary  condition, 
such  as  X  =  0  or  F  =  0,  would  be  known  at  each  end  of  the  system,  then  the  equations  can  be  solved  to 
yield  the  natural  frequencies  and  principal  modes  of  vibration.  The  precise  manner  in  which  this  would 
be  carried  out  is  not  readily  apparent,  and  wiU  be  iUustrated  below.  Two  methods  are  suitable.  The  first 
yields  the  nth-order  characteristic  equation,  which  can  then  be  solved  for  the  n  natural  frequencies.  Each 
natural  frequency  can  then  be  substituted  into  individual  stages  of  the  determinations  to  give  the  amplitude 
ratios  and  thus  the  modal  pattern.  This  procedure  is  feasible  only  if  the  number  of  degrees  of  freedom 
is  small. 
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The  second  method  is  to  follow  a  numerical  procedure  such  as  Holzer’s  method.  By  assuming  a  value 
for  (0^  and  assigning  a  value  to  the  unknown  boundary  condition  at  station  0,  then  the  unknown 
parameters  of  x  and  F  can  be  determined  at  Station  1  by  the  related  matrix  multiplication.  This  process 
can  be  continued  from  one  station  to  the  next  until  the  fmal  station  n  is  reached  and  the  boundary 
condition  there  is  checked,  yielding  the  error.  A  new  value  of  aP'  would  then  be  assumed  and  the  entire 
procedure  would  be  repeated,  resulting  in  a  new  error  value.  The  process  can  be  continued  until  the  error 
is  brought  to  zero — or  rather,  to  acceptable  small  value. 

To  illustrate  the  above  methods,  consider  the  mass-spring  system  shown  in  the  Figure  5.  Let  us  now 
use  transfer  matrices  to  detennine  the  principal  modes  and  natural  frequencies  by  the  first  method 
described  above.  Assume  kj  =  kj  =  k3  =  k  and  mj  =  m,  mj  =  2m. 


(0)  (1)  (2)  (3) 


Figure  5.  Mass-soring  system  for  transfer  matrices. 

The  boundary  conditions  for  Sta.  0*^  are  X  =  0,  F  =  F^,  and  for  Sta.  3*^  are  x  =  0,  F  =  F3,  where  F3 
is  unknown.  The  first  transfer  matrix  would  be  written  for  Sta.  1*^  in  terms  of  Sta.  0*^  as 


a.) 


1 


-m©^ 


1 


1 


k 


m©^ 


(172) 


(173) 
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Next,  the  transfer  matrix  for  Sta.  2  is  written  as  follows: 
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The  computations  in  the  first  column  of  certain  of  the  transfer  matrices  shown  by  the  vertical  line  are 
omitted  as  they  are  not  needed,  due  to  the  zero  in  the  first  row  of  the  state  vector  for  Sta.  0*^.  Continuing, 


From  the  first  algebraic  equation  of  this  final  matrix,  we  have 


(175) 


(176) 
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o 


hence 


mV 


(177) 


which  represents  the  characteristic  equation  and  has  the  roots 


mV  ^  3  -  vT  3  +  v/s" 


(178) 


This  defines  the  natural  frequencies  as 


2  _  3  -  v/T 


®1  = 


m 


2  3  +  v/3  k 


(O2  = 


m 


(179) 


Setting  Xj  =  1  and  also  substituting  ©j  into  Equation  (173)  gives 


Fo*k.  F, 


V^-1 


k. 


(180) 


Then  substituting  these  into  Equation  (174)  gives 


■  _  1  + 
2 


(181) 


Similarly,  setting  Xj  —  1  and  substituting  0)2  into  Equation  (173)  gives 


F„  =  k,  F,  .  Zl/j-  ■  k 


(182) 
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and  substituting  these  into  Equaticm  (174)  yields 


X2  = 


1-V^ 

2 


Thus  the  principal  modes  are  as  follows: 


2  _  3  -  k  _ 


COj  = 


,, _ =  0.6334  A,  *  1,  A,  _ 

2  m  m  2 


=  ^  =  1. 


366 


(183) 


2^3  ^  2.366  JL,  Aj  =  1,  Aj  =  ^  ~  ^  =  - 

2  m  m  ^2 


“2  = 


0.361 . 


(184) 


For  a  torsional  system  composed  of  disks  on  a  massless  shaft,  the  analysis  and  formulation  of  the  transfer 
matrices  are  identical  to  that  of  the  rectilinear  mass-spring  system.  Thus  considering  the  torsional  systems 
and  corresponding  fiee-body  diagrams  shown  in  Figure  6,  the  transfer-matrix  relations  for  the  torsional 
shaft  (spring)  Kj  are 


i:- 


1  I 

k 

0  1 


-ij 


e  1R 

M 


‘Jj-l 


(185) 


or  in  general  form 


(186) 


For  disk  the  matrix  relations  for  harmonic  motion  ate 


f  e  1  *  f  1  o]  /  e  1  >- 
r-Jj  iJjHlj’ 


(187) 


having  the  general  form 


(188) 
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Figure  6.  Torsional  systems  associated  with  gear  box  of  Ml  turbine  engine. 


Substituting  Equation  (186)  into  Equation  (188)  gives 


(189) 


where 


1 

K 


1 

V 


10)2 


(190) 


If  Equation  (189)  is  applied  to  successive  stations  of  the  torsional  system,  then 


=  [Q]„,  W 


(191) 


36 


where 


[Q]„!  =  [P]„[F]„[P]„_  i[F]„_  1  ...  [P]i[F]j.  (192) 

Equations  (185)  through  (192)  are  the  same  as  those  for  the  mass-spring  system  except  that  m  has  been 
replaced  by  I,  X  has  been  replaced  by  0,  and  the  k  now  represents  a  torsional  spring  constant 

Let  us  now  consider  a  torsional  system  shown  in  Figure  7  and  let  us  employ  transfer  matrices  to 
determine  the  principal  modes  and  natural  frequencies  by  the  second  method  described  above  (the  Holzer 
procedure). 


Figure  7.  Torsional  system  of  Ml  gas  turbine  engine  with  bearing. 


The  boundary  conditions  for  Sta.  0*^  are  6  =  0,  M,  =  M^,  =  1  and  for  Sta.  3*^  arc  0  =  0,  Mj  =  0,  where 
03  is  unknown.  Assuming  the  first  trial  frequency  value  as  co^  =  1.0  (K  I),  and  referring  to  Equation  202, 
the  first  transfer  matrix  for  Sta.  1^  in  terms  of  Sta.  0^  would  be  written  as 


a.) 


1 

2k 

1 


0 

t 


(193) 
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As  in  the  preceding  example,  the  elements  of  the  first  row  of  the  transfer  matrix  shown  by  the  vertical 
line  are  not  included  as  they  are  not  needed,  due  to  the  zero  in  the  first  row  of  the  state  vector  for  Sta.  0^. 
This  condition  also  applies  to  subsequent  transfer  matrices  in  this  example. 

Next,  the  transfer  matrix  for  Sta.  2*^  is  written, 

b.) 


38 


Next 


c.) 


1 

6k 

_£ 

"6 


R 

0 


lx±.Ix 
6k  k 


-k  X  —  +  0  X 
6k 


R 

0 


(195) 


From  the  second  algebraic  equation  of  this  matrix,  Mjj  =  -2  representing  the  error  since  M.  should 

6  ‘ 

be  zero. 

A  new  value  of  is  now  chosen  and  the  process  is  repeated,  resulting  in  a  new  error.  This  is 
continued  until  the  Mj3  error  is  reduced  to  an  acceptable  value,  and  the  corresponding  eo^  is  then  a  correct 

natural  frequency  value.  Plotting  the  error  vs.  the  value  aids  in  interpolating  for  the  correct  frequency 
values.  Such  continued  calculations  are  not  shown  here,  but  tire  correct  (acceptable)  aP-  values  can  be 
verified  as 


©5  =  0.2168,^,  ©2  =  1.0965  2,  ©,  =  2.1035  2. 

I  ^  I  ^  I 


(196) 
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The  corresponding  errors  are 


-  0.00000,  =  0.0013,  -  0.00023. 


(197) 


The  principal-mode  amplitude  ratios  can  be  determined  by  substituting  the  natural  frequency  values  into 
the  transfer  matrix  equations  for  each  station. 
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APPENDIX: 

TYPICAL  SINGLE-DEGREE-OF-FREEDOM  DIFFERENTIAL  EQUATIONS 

AND  THEIR  SOLUTIONS 
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1.  FREELY  VIBRATING— NO  DAMPING 


Solution 


where 


X  +  (O^X  =  0 


X  =  Xf,  cos  ©t  = _ sin  ©t 

© 


Xq  -  initial  displacement  -  in 

Vq  -  initial  velocity  -  in/s 

©  -  natural  circular  fiiequency  -  rad/s 

2.  FREELY  VIBRATING— WITH  VISCOUS  DAMPING 


Solution 


X  +  2^©X  +  ©^X  =  0 


X  »  e  Ajcos(©o  ^  “  4*) 
=  ©  [l  -  ^  ^ 


© 


Aj  = 


Xq  + 

^0  +  Xo^©^ 

2 

0 

o 

3 

1 

T 


(j)  =  tan 


-1 


Vq  +  ^©Xq 


^  Xo©o  ^ 


Xq  -  initial  displacement  -  in 
Vq  -  intial  velocity  -  in/s 

©Q  -  damped  natural  circular  frequencies  -  rad/s. 
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3.  CONSTANT  FORCED  VIBRATION~NO  DAMPING 


X  +  co^x  =  G 


X  =  Xq  cos  tot  + 


^0  .  ,  G 

_  smcot  + _ 


to 


to 


where  G  =  a  constant  independent  of  time. 


4.  FORCED  VIBRATION  WITH  NO  DAMPING 


X  +  to^X  =  G(t)  =  GflCOsOt 


Solution 


X  =  XoCOStOt 


Vo  . 

+  —  sintot  + 
to 


to^  - 


[cosOt  -  costot 


O  -  forcing  frequency  -  rad/s 


Steady-state  solution 


X 


s 


to^  - 


cosOt 


Gq  1 
1  - 


cosQt , 


where 


r  = 


Q 


to 
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Subcases  (steadv-state  solution) 


(1)  Go  =  Gia)2 

Gi 

X,  =  _ cosAt 

*  l-r2 


(2)  Go  =  G2ft^ 

r2 

X3  =  G2 _ cosflt 

1  - 


5.  FORCED  VIBRATION  WITH  DAMPING 


X  +  2^(dX  +  (O^X  =  G(t)  =  GoCOsOt 


Solution 


X  =  Aje^“‘  cos  |(Bot  - 


(^t  -  <j)jJ 


(1  -  rf  *  (25rP, 


1 

T 


+ 


+  ^coXo 
Xo“o 


) 


1 

T 


<l> 


=  tan"^ 


f^o 


+  ^caXo' 

XqCOq 


<|)i  =  tan"^  . 

V 


25r  1 
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Steady-state  solution 


Xs  = 


cos 


(^t  -  $l) 


Subcases  (steadv-state  solutions') 


(1)  Go  =  Gi(o2 


X,  =  . 

Gicos^flt  -  (t)j) 

s 

[(l  -  rf  +  (2^0^]^ 

(2)  Go  =  GjQ^ 

X  =  _ 

G2T^  cos^Qt  -  <|)i) 

S 

(l  -  rf  .  (2|r)^j^ 
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GLOSSARY  -  NOMENCLATURE  AND  UNITS 


k  -  spring  constant  or  spring  modulus  (linear-lb^n)  (toisional-in/lb/rad) 
m  -  mass-lbrin/s^ 
mg  -  w-weight-lb 

J  -  polar  mass  moment  of  inertia-lb/in/s 
I  -  area  moment  of  inertia-in^ 
g  -  gravitational  constant-386.4  in/s^ 

C  -  viscous  damping  constant-(linear-i^^)  (torsional-  ) 

in  rad 

c  -  damping  ratio  =  C/Cg 

-  critical  viscous  damping  constant — that  amount  of  damping,  which  if  introduced  into  the  system  just 
prevents  the  system  horn  vibrating 

X  -  displacement-in 
X  -  velocity-in/s 

X  -  acceleration-  in/s^ 

6  -  angular  displacement-radian 

6  -  angular  velocity-rad/s 

6  -  angular  acceleration-rad/s^ 
t  -  time-seconds 
CO  -  natural  circular  Itequency-rad/s 
Q  -  forced  circular  frequency-rad/s 

r  =  ^  ratio  of  forced  to  natural  circular  frequency 
CO 
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f 


0) 

2n 


frequency  of  vibration  in  cycles/s 


5  =  - =  logarithmic  decrement  =  In 

[i  - 

L  =  Polar  area  moment  of  inertia-in'* 
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